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Abstract 

When a diquark does not fragment directly but breaks in such a way that only one of its 
quarks gets into the produced baryon, the latter is produced closer to mid rapidities. The 
relative size of this diquark breaking component increases quite fast with increasing energy. 
We show that at a given energy it also increases with the atomic mass number and with the 
centrality of the collision and that it allows to explain the rapidity distribution of the net 
baryon number (p-p) in SS central collisions. Predictions for Pb-Pb collisions are presented. 
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In most string models particle production in a pp collision occurs mainly in the form 
of two diquark-quark strings. There is strong experimental evidence [|lj that with a large 
probability the diquark acts as a single entity and fragments directly into a leading baryon. 
Since the diquark is fast in average, taking a large fraction of the incoming proton momen- 
tum, the produced baryon will be in the proton fragmentation region. Moreover, it is easy 
to see that the corresponding baryon spectrum at y* ~ decreases, quite fast with energy 
- roughly as s^ 1 H-[0|- A diquark breaking component in which only one of two quarks 
gets into the produced baryon is of course present. This component produces baryons at 
smaller values of y* and the produced baryon distribution at y* ~ only decreases like 
s- 1 ^ 0. In the Dual Parton Model (DPM) and Quark Gluon String Model (QGSM) such 
a component is not explicitly taken into account ||, [|J. However, it has been shown by 
Kopeliovich and Zakharov || that the energy dependence of da/dy at y* = for the dif- 
ference p minus p measured by two different collaborations 0, in the ISR energy range 
can only be understood with a diquark breaking component of a relative size of about 20 % 
which has to be introduced in the string models as a separate component. The authors of 
ref. have also proposed an interpretation of the diquark breaking component as a result 
of a color exchange which converts the diquark from a {3} to a {6} color state. While this 
is an attractive mechanism, it may not be the only one. If to assume that the relative longi- 
tudinal momenta of the quarks inside the diquark may be substantially different, a diquark 
breaking with excitation of a high-mass state could also result in the baryon production at 
mid rapidities. 

In the present work we do not specify the origin of the diquark breaking component 
and do not attempt to compute its size theoretically (for a calculation in perturbative QCD 
see ref. 0). Following the phenomenological analysis in ref. || we determine the latter 
from experiment. Also, we concentrate ourselves to present CERN energies of \/s ~ 20 
GeV and make only some comments about our expectations for RHIC or LHC energies. We 
shall show that in heavy ion collisions the relative size of the diquark breaking component 
increases with atomic mass number and with the centrality of the collision. We show that 
in central SS collision the rapidity distribution of the difference p-p - which in DPM has a 
sharp dip at mid rapidities || not present in the data [0] - is consistent with experiment 
after introducing the diquark breaking component - with relative size determined from pp 
data as explained above. The corresponding distributions for central Pb-Pb collisions at 
present CERN energies can then be computed, and are given as a prediction of our diquark 
breaking mechanism. 
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In DPM, the dominant component in pp scattering is the two string configuration de- 
picted in fig. 1. Here the diquark fragments directly into a final state baryon. The rapidity 
distribution dN / dy of a hadron in a string is obtained from a convolution of the momentum 
distribution and fragmentation functions PJ. All formulas and parameters relative to baryon 
production in nucleon-nucleon collision are given in refs. and ||. The corresponding re- 
sults will be denoted dNpp/dy (where DP stands for diquark preserving). As mentioned 
above there is also the possibility that the diquark breaks in the way depicted in fig. 2. 
In this case the baryon rapidity distribution follows that of the valence quark and thus the 
baryon is slower than in the DP component of fig. 1. We then have 




Fig. 1 Fig. 2 

Figure 1: Diquark preserving component in nucleon-nucleon scattering. 
Figure 2: Diquark breaking component in nucleon-nucleon scattering. 



~^y~^ = J \?Qv(y) + Piv{-y)\ > (!) 
where DB stands for diquark breaking and p qv (y) is the valence quark rapidity distribution. 
This function and the value of the constant C are given below. Obviously, in the final 
formula, the DB component has to be weighted by the ratio erDs/oin of the DB over the 
inelastic cross-section. 

Let us now turn to a nucleon-nucleus (NA) collision. Here we can also have a DP and 
DB components. In the case of a double inelastic collision they are depicted in figs. 3 and 
4, respectively. In the NA case (as well as in a NN configuration involving at least four 
strings) it is also possible to have a configuration as the one depicted in fig. 5 in which 
the baryon number is associated to a gluon or sea quark of the incoming nucleon. For the 
moment we disregard such a configuration but shall come back to it later on. Let us now split 
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Fig. 3 Fig. 4 

Figure 3: Diquark preserving component for nucleon-nucleus scattering 
for two inelastic collisions. 

Figure 4: Same as fig. 3 for the diquark breaking component. 



the nucleon-nucleon (NN) cross-section into its diquark preserving and diquark breaking 
components a in = a DP + a DB . We assume that once the diquark has been destroyed in a 
collision with one nucleon of the nucleus it cannot be reconstructed in further collisions with 
other nucleons of the nucleus. The NA cross-section involving n inelastic NN collisions is 
then given by 

<tn(b) = (^j t (") ^ Tlib) [1 - o m T A (b)) A ~ n . (2) 




Figure 5: Same as fig. 4 but here the produced baryon has the rapidity 
distribution of a sea quark. 
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Here is the standard nuclear profile function at impact parameter b, normalized to 

unity. In (2) we have replaced the usual factor a™ n T\ corresponding to the cross-section for 

n inelastic collisions by the product ( 1 f) cr DB a liyp T/L = ( a ? n ~ °"Sp)^a(^)- Indeed only 

i=i 

the term a^ P T\ will contribute to the diquark preserving cross-section. Summing in n we 
have 

«) = EO) = l-[l-^T A (6)f . (3) 

71=1 

Eq. (3) shows that the diquark preserving cross-section belongs to a class of process || 
which has only self-absorption (or self-shadowing). Obviously 

= °Z A Q>) ~ ■ (4) 

Since <Jdb < &im h is clear from (3) and (4) that a pip increases faster with A than &pp. 
Actually when a flB is sufficiently small to neglect in (3) second and higher powers of <Jdb, 
a DB wm increase linearly with A. This proves the result stated above that the relative size 
of the DB component increases with increasing A. The result can be easily generalized to 
an AB collision. We have P| 

a* B B {b) = l-{l-a DB T AB {b)) AB . (5) 

where Tab(P) = J d 2 s Ta(s) T B {b — s). For a DB sufficiently small we get from (5), after 
integration in impact parameter, a^p = AB ap>B- 

In order to compute the nucleon rapidity distribution we only need to specify the values 
of (Xob as well as the value Cnn of the constant C in (1) for an iviV collision, as well as 
the function p? v (y) = Z Pg v (Z) in (1). Here Z = e~ Ay , Ay = y — Ymax- The function 
Pq v {Z) is explicitly given in DPM and QGSM and depends on the number n of inelastic 
collisions suffered by the nucleon {n = 1 in Figs. 1 and 2 and n = 2 in Figs. 3-5). Indeed, 
in these models the nucleon splits into 2n components (valence quark, diquark and n — 1 
gluons or q s -q s pairs) and the joint moment distribution is written as a product of the 2n 

In 

factors controlling the Z ~ behaviour of the 2n constituents times 5(1— J2 Zi)- In QGSM, 

i=i 

where both valence and sea quarks have (at the energies considered here), the same Z -1 / 2 
behaviour, one has 

12 12 12 I ^ \ 

p n (Zi, ■ ■ ■ Z 2n ) = C n Z l Z 2 ■ ■ ■ Z 2n _ x Z 2n 5 |1 — Zij 
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where Z 2n corresponds to the diquark and the constant C n is determined from the norma- 
lization to unity. In this case one gets by integrating over all variables except the valence 
quark one 

p» q M = z pi(z) = c n z^(i-zr + ^ . (6) 

In DPM, on the other hand, one assumes a Z~ l behaviour for the sea quark with Z = 



Z 2 + p 2 / s. This only changes the effective power of (1 — Z) in (6) and we have checked 
that at present energies the differences are only of a few per cent. 

The quark momentum distribution in (6) is also appropriate for the diquark breaking 
case when the latter takes place via the color exchange mechanism described in ref. ||. 
Indeed, in this case the diquark, while changing its color, retains essentially its momentum. 
However, in the case of diquark breaking via longitudinal momentum transfer, where the 
two quarks in the diquark are separated, the number of independent constituents is now 
2n + 1 (see figs. 2 and 4 for the cases n = 1 and n = 2 respectively). Moreover, the X- value 
of one of the quarks has to be larger than Z - i.e. the momentum fraction of the quark 
carrying the baryon number. 
One then has instead of (6) 



, r 1 - 2 dX 1 
ftM^Z pl(Z)=C> n Z^j z —(i-X-Z)^ . (7) 



'i-z dX 
lz ~JX 

In the numerical calculations we have used eq. (7). We have checked that using (6) the 
changes in our results are at most of the order of 10 %. 

For the constant C in eq. (1) we take, in the case of a NN collision, C = 0.84. The 
factor 0.84 originates as follows. In NN collision we produce a net baryon number equal 
to 2 with equal number of protons and neutrons. Moreover, at the present energies we 
know that the number of A/S° is about 0.1 while the number is about 0.06 0. 



For the diquark splitting component of figs. 2 and 4, however, the production of strange 
baryons is twice as large as in the case of direct diquark fragmentation since the s-quark 
can be created in either side of the valence quark. (In the component depicted in fig. 5 the 
production of strange baryons would be enhanced by a factor 3). Therefore from baryon 
number conservation the average multiplicity of either protons and neutrons will be 1 - 0.1 
- 0.06 = 0.84 (remember that p^ v (y) is normalized to unity). 
Finally we take Cj n = 32 mb and from ref. 0^ = 7 mbQ. 

■''Actually, two values of (TdBi dependent on the form of the proton wave function, are quoted in 0]. 
The lower one, describes well the experimental data on the energy dependence of (dN/dy) p — (dN/dy) p 
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We are now ready to compute the nucleoli rapidity distribution in NN and AB collisions. 
We concentrate on the difference Ap = p — p which has been measured by the NA35 
collaboration in peripheral and central SS collisions ]7). It is given by 

iv) = (l-—)^ + (y) . (8) 



dy V a in ) dy a in dy 

Here dN DB /dy is given by eqs. (1) and (7), with n = 1, while dN DP /dy has been computed 
in ref. (all formulae and numerical constants are the ones given there). 
Likewise in an AA collision we have 



dN AA^Ap 

.(y) = (n A -n D A B ) 



7 i\tNN — >Ap 

dy 



+ n°A B 



n=n/nA 



dy 



• (9) 

n=n/nA 



dy 
where 

n D A B {b) = A (1 - a DB T AA (b)) A / a^(b) 

and n A has the same expression except that a^B is replaced by <Tj n - both in the numerator 
and in the denominator. Eq. (9) is easily generalized for the case of collision of different 
nuclei A and B. Distribution function dNns/dy in (9) is given by eqs. (1) and (7) with 
n = n/n A , i.e. the average number of collisions per participant nucleon. Again dN^pjdy is 
the one computed in ref. . For central collisions, n A and n A B in (9) have to be computed 
for b ~ 0. 

Note that by construction eqs. (8) and (9) reproduce the rapidity distributions of ref. [§J 
for odb — 0. The extra nucleon stopping for a^B 7^ is due to the fact that n A B increases 
much faster than n A with increasing A - and to the different shapes of DP and DB rapidity 
distributions. The former has a pronounced dip at y* ~ 0, while the latter is concentrated 
at mid rapidities. 

The results for peripheral and central SS collisions at CERN energies are shown in fig. 
6 and compared with the data from the experiment NA35 0. We also show the results 
obtained with o DB = ||. We see that the dramatic dip present in the latter case for 
central SS collision has been largely filled in by the contribution of the DB component in 
agreement with experiment. The prediction for a central Pb-Pb collision is also shown. In 
this case the dip is converted into a broad plateau. 



at y* ~ 0, (line II in fig. 3 of that reference). This line is obtained from our eq. (1), when weighted by 
odb I <^m with (7db ^6-^7 mb. In the notations of ref. Q| this line corresponds to <jdb = 3.5 mb/K with 
K ~ 0.5. In the following we shall use the value odb = 7 mb. (Using odb = 6 mb the changes in our 
results are less than 5 %). 
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We conclude with two remarks on strange baryon production and on the extrapolation 
of our results to much higher energies. Concerning the first, we have already explained that 
the ratio of strange over non-strange baryons is twice as large in the DB component as in 
the DP one : the corresponding A enhancement is, however, rather small. We obtain an 




Figure 6: Rapidity distribution of the p-p difference in AA collisions 
at y/s = 20 GeV. The thick dashed curve corresponds to the diquark 
preserving mechanism of ref. [4] for central SS collision. The thick solid 
curve shows the result of our calculations including the diquark breaking 
component for central SS collision. The full circles are the data points 
from [7]. 

The same for peripheral SS collisions is shown by the thin solid and 
dashed curves and the square points. Here the theoretical curves are for 
AW collisions, normalized to the data. 

The dotted curve shows our prediction for central Pb-Pb collision scaled 
down by a factor 1/5. It corresponds to nucleon minus antinucleon - 
rather than to proton minus antiproton. 
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extra A average multiplicity of 1.1 in central SS and 13 in central Pb-Pb (with extra y* — 
densities of 0.3 and 3.8, respectively). As discussed in ref. [|J the existing data require final 
state interactions such as it + N —>■ K + A. The results obtained with the nucleon rapidity 
distributions computed above are close to the ones obtained in ref. [|J, where some extra 
stopping had been arbitrarily introduced (see full lines in fig. 3 of ref. Q]). More precisely 
the A distribution in central SS is increased by about 10 % as compared to the one given 
in M, and the one for central Pb-Pb collisions is increased by about 15 % at y* ~ but is 
narrower than the one given in 0]. 

Finally, we turn to the DB component in fig. 5 which has not been considered so far. As 
emphasized in ref. in this component the baryon number follows the distribution of a sea 
quark or gluon and would lead to a baryon distribution flat in rapidity and independent of 
energy at very high energies. At present CERN energies, however, such a component would 
produce essentially no change in our results. Indeed, as discussed above, a sea quark has an 
Z^ 1 distribution in DPM, and we have checked numerically that at the energies considered 
here this distribution differs by less than 10 % from the one of a valence quark. However, 
at energies of RHIC and especially LHC this component would provide an interesting and 
efficient mechanism of nucleon stopping. As discussed above it would also be efficient in 
producing strange baryons at mid-rapidities. Actually two extreme scenarios may be con- 
sidered : 1) the component depicted in fig. 5 is strongly suppressed for some dynamical 
reason and can be neglected and 2) at RHIC energies and beyond, this component is the 
dominant part of DB - and is independent of s. In both cases the DP component at y* ~ 0, 
which falls roughly as s -1 , is already negligible at RHIC energies. As for the DB one, it 
falls in the first scenario as s -1 / 4 giving a y* ~ net baryon number, p-p, of 10 at RHIC 
and 2 at LHC. In the second scenario, one gets a net baryon number at y* ~ of 10 at both 
RHIC and LHC, which is about 30 % of the one predicted at 160 GeV/c (see fig. 6). 

Finalizing this work we learned that D. Kharzeev has just published paper fll| on a 



similar subject, and we would like to compare our results. The very observation of [|TTJ that 
gluons can trace baryon number was already done in ||, and in addition the cross section 
was evaluated. We disagree with JTT[| , however, on the rapidity dependence for the gluonic 
mechanism of baryon number flow. Solid arguments presented in [|1^, |T^] (see also ||14j| ) 
and in the present paper prove that this mechanism is rapidity independent. However, this 
point is not important for phenomenology at present energies since the data [||, [7J are 
dominated by the preasymptotic mechanism of a valence quark exchange accompanied with 
the diquark destruction (exchange of a valence quark together with the string junction in 
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terms of topological expansion [T5|]). Finally, in contrast to fTI |, we have presented a full 
calculation of baryon stopping in AA collisions. 
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